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We study the excitation of axial quasi-normal modes of deformed non-rotating black holes by 
test-particles and we compare the associated gravitational wave signal with that expected in general 
relativity from a Schwarzschild black hole. Deviations from standard predictions are quantified by 
an effective deformation parameter, which takes into account deviations from both the Schwarzschild 
metric and the Einstein equations. We show that, at least in the case of non-rotating black holes, it 
is possible to test the metric around the compact object, in the sense that the measurement of the 
gravitational wave spectrum can constrain possible deviations from the Schwarzschild solution. 


I. INTRODUCTION 

So far, general relativity has successfully passed all 
the experimental tests Jj. However, its most interest¬ 
ing predictions are still to be verified. The theory has 
been tested in weak gravitational fields, mainly with pre¬ 
cise experiments in the Solar System and accurate radio 
observations of binary pulsars. The agreement between 
theoretical predictions and observational data is today 
confirmed with a precision ranging from a few percent 
to about 10 -5 . There is now an increasing interest to 
test the theory in other regimes, in particular at very 
large scales and in strong gravitational fields. Tests of 
general relativity at very large scales are relevant for the 
problems of dark energy and dark matter: in the last 20 
years, this has been a very active research field [2]. The 
ideal laboratory to test strong gravity is the spacetime 
around astrophysical black hole (BH) candidates (H [4] . 
The spacetime geometry around these objects has still 
to be verified and there is a number of theoretical argu¬ 
ments suggesting that macroscopic deviations from stan¬ 
dard predictions are possible m- 

In 4-dinrensional general relativity, uncharged BHs are 
described by the Kerr solution and they are completely 
specified by only two parameters, namely the mass M 
and the spin angular momentum J. This is the well- 
known no-hair theorem El- Astrophysical BHs should 
form from the collapse of very massive stars, after the 
latter have exhausted all their nuclear fuel m ■ Accord¬ 
ing to the theory of general relativity, the final prod¬ 
uct of the collapse should be well described by the Kerr 
solution. Initial deviations from the Kerr metric are 
quickly radiated away through the emission of gravita¬ 
tional waves m ■ For macroscopic objects, the equilib¬ 
rium electric charge is extremely small and is reached 
soon because of the highly ionized host environment m- 
The effect of an accretion disk is normally negligible, as 
the disk mass is many orders of magnitude smaller than 
the mass of the central object (13] . 
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Compact objects in X-ray binaries are classified as 
stellar-mass BH candidates if they exceed 3 M 0 , which 
is the maximum mass expected for neutron and quark 
stars M- The supermassive objects at the center of 
every normal galaxy are called supermassive BH candi¬ 
dates because they cannot be explained as clusters of 
non-luminous objects: the expected cluster lifetime due 
to evaporation and physical collisions would be shorter 
than the age of these systems US). The non-observation 
of thermal radiation emitted by the possible surface of 
these objects may also be interpreted as an evidence for 
an event/apparent horizon |I^ (but see Refs. JT, TBj and 
there are also scenarios in which no thermal component 
should be expected 5, 6 ). In the end, both stellar-mass 
and supermassive BH candidates can be naturally inter¬ 
preted as the Kerr BHs of general relativity and they 
could be something else only in presence of new physics. 
However, a direct observation confirmation is completely 
missing. 

The nature of BH candidates can be potentially tested 
by studying the properties of the radiation emitted by 
the accreting gas [19l - t26] . The electromagnetic spectrum 
of these objects has indeed features determined by the 
motion of the gas in the accretion disk and by the prop¬ 
agation of the photons from the disk to the distant ob¬ 
server. The study of these features can potentially con¬ 
strain the spacetime geometry and check whether these 
objects are the Kerr BHs of general relativity. Current 
observations cannot do it, mainly because there is a fun¬ 
damental degeneracy among the parameters of these sys¬ 
tems: the relativistic features produced in a non-Kerr 
background cannot be distinguished from those produced 
around a Kerr BH with a different spin [??H3(J] . With the 
available data, we can rule out some exotic BH alterna¬ 
tives, like some kinds of compact objects ins eh and 
of traversable wormholes [33]. Non-Kerr objects with a 
horizon are much more difficult to test. Current con¬ 
straints on stellar-mass BH candidates from the disk’s 
thermal spectrum are reported in [33, 35) ■ The possi¬ 
bility of testing the Kerr metric in the future seems to 
depend on the possibility of combining several observa¬ 
tions of the same source to break the degeneracy among 
the parameters of the system El [37]. 
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A complementary approach, not yet available, is to use 
the gravitational wave signal [38H40] . This approach has 
advantages and disadvantages. Assuming geodesic mo¬ 
tion, the properties of the radiation emitted by the gas in 
the accretion disk only depend on the background metric, 
and therefore they can be used to test the Kerr metric. 
They cannot distinguish a Kerr BH in general relativity 
from a Kerr BH in an alternative theory of gravity SB- 
On the contrary, the emission of gravitational waves de¬ 
pends on both the background metric and the field equa¬ 
tions of the theory, with the advantage that it may be 
possible to distinguish Kerr BHs in different theories and 
the disadvantages that a more rigorous treatment would 
require the knowledge of the field equations. However, 
a phenomenological description of the problem is possi¬ 
ble and it has been already employed, for instance, in 

Refs, mm- 

In this work, we investigate the gravitational wave sig¬ 
nal emitted when the axial quasi-normal modes of a BH 
are excited by the passage of a test-particle. As an ex¬ 
ploratory work, we consider the simplest case of non¬ 
rotating BHs. We derive the master equations and we 
compute the wave form and the energy spectrum as a 
function of the deviations from general relativity and 
of the energy and the angular momentum of the test- 
particle. Deviations from standard predictions are quan¬ 
tified by an effective deformation parameter, which is 
used to take into account a possible non-Schwarzschild 
background and corrections to the Einstein equations. 
We find that the gravitational wave signal is mainly 
determined by the geodesic motion of the test-particle, 
while the contribution from the proper excitations of the 
spacetime is smaller. 

The content of the paper is as follows. In Section [TTJ 
we review the excitation of axial quasi-normal modes of a 
Schwarzschild BH by a test-particle in general relativity. 


In Section |III[ we consider the same phenomenon in the 
case of a non-rotating deformed BH. Since we compute 
the gravitational wave signal from the Einstein equations, 
eventually our predictions are parametrized by an “ef¬ 
fective” deformation parameter. Such a deformation pa¬ 
rameter is not the same as the one appearing in the back¬ 
ground metric and therefore it cannot be directly com¬ 
pared with the deformation parameter constrained by the 
studies of the electromagnetic spectrum. In Section [TV} 
we present the results of our calculation and we com¬ 
pare the gravitational wave signals from Schwarzschild 
and deformed non-rotating BHs. Summary and conclu¬ 
sions are reported in Section |V| Throughout the paper, 
we use units in which Gn = c = 1. 


II. SCATTERING OF PARTICLES BY A 
SCHWARZSCHILD BLACK HOLE 

To study the gravitational wave signal by a test- 
particle perturbing the Schwarzschild background, we 
can proceed as follows [44]. At the zeroth order, the 
Einstein equations are R° a g = 0 and the unperturbed so¬ 
lution is the static Schwarzschild background g^g- The 
first order term in the perturbed metric is g] y g and it can 
be derived by 

RU = (r aP - . ( 1 ) 

We expand the metric g\g in tensor spherical harmonics. 
In this work, we only consider axial perturbations. The 
equations governing polar perturbations can be obtained 
from the axial ones by a coordinate transformation. We 
choose the so-called Regge-Wheeler gauge and the per¬ 
turbed line element reads [44] 


ds 2 = gl f} dx a dx 13 + 2 sin6> y ^ lm ^ ^ [/ Zm (f, r)dtd(t> + h lm {t , r)drd(f)} 


l,m 


r)dtdd + him(t, r)drdO} . 


sin 9 ' 86 

l,m 


( 2 ) 


The matter source in the Einstein equations is given by the stress-energy tensor of a point-like particle moving along 
the geodesics of the Schwarzschild background, namely 

dT d7 a d~P 1 

T< = ^~df id ^ (r " m) 6 {e - ( 3 ) 

where g is the particle’s rest-mass, z a = (T(f), f?(i), 0(f), <I>(f)) is the world-line trajectory of the particle in term of 
the Schwarzschild t-coordinate, and r is the particle’s proper time. The tensor T a d is expanded in tensor spherical 
harmonics as well [45] . After some tedious calculations, the components t6 , r</>, and 96 of Eq. (|T|) provide the following 
equations 


r ) - + Vi(r)Z lm (t,r) = Si m (t,r ), 

d d 

= — [r*Z im (f,r)] , 


(4) 


( 5 ) 
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where 


Z,„(f,r)=( 


and r* is the tortoise coordinate, namely 


r* = r + 2M In 


(~2M 


-1 ■ 


The potential Vi(r) and the source term Si m (t,r) are given, respectively, by 


y I (r)=(l- — 


1(1 + 1) 6 M 

r 2 r 3 


where 


Sim(t,r) = -167 ri ( 1 - —) [r® 


Dim — 
Qlm — 


fim 


r J ( dr 


Li l 


, 2 M\ „ 

1 -— J Dim 


0 «- 


l(l + l)(l-l)(l + 2) E r 4 V r / 
i[i 


- S(r- R(t )) 


50 


o=f,$(t) 


Z(Z + 1 ) r 2 (r — 2 M) dr 


A* 

<5 (r - U(r)) 


<90 


©=i>*W 


( 6 ) 

(7) 

( 8 ) 

(9) 


( 10 ) 

( 11 ) 


Here 1? and L z are, respectively, the specific energy and the axial component of the specific angular momentum of 
the particle, which is assumed to move on the equatorial plane 0 = ^. The wave form of the signal is Zi m (t) = 
limr^+oo Zi m (t,r). The energy spectrum of gravitational waves at infinity is proportional to 


dW 

dui 


+00 l 

m + 1){l + 2) \ z ^)\ 2 , 

1—2 m=—l 


( 12 ) 


where Zi m (uj) is the Fourier transform of the wave form 
Zi m (t). We remind the reader that there are no monopole 
(l = 0 ) and dipole (l = 1 ) contributions in general rela¬ 
tivity, and that the leading order term is the quadrupole 
(1 = 2) one. Moreover, multipoles with m = l — 1, 1 — 3, 
etc. correspond to axial perturbations, while those with 
m = l, l — 2 , etc. correspond to polar perturbations. As 
a preliminary study, here we are restricting our attention 
to axial perturbations only, and therefore we will ignore 
the signal associated to polar perturbations. 


III. SCATTERING OF PARTICLES BY A 
DEFORMED BLACK HOLE 


In the case of an alternative theory of gravity, one could 
proceed in the same way, namely starting from the corre¬ 
sponding BH solution and considering the perturbations 
generated by a test-particle moving along the geodesics 
of that metric (in the case of a metric theory of gravity). 
The evolution of these perturbations and the associated 
gravitational wave signal at infinity are then determined 
by the field equations of the theory. The final result 
depends on both the background metric and the field 
equations. In this sense, the approach is more powerful 
than the study of the properties of the radiation emit¬ 
ted by the gas in the accretion disk, as the latter is only 


sensitive to the background metric. However, if we use 
gravitational waves, we can only test a particular theory, 
while the actual fundamental theory may not be known. 

In this paper, we explore a different approach. We con¬ 
sider the Johannsen-Psaltis metric [32], which describes 
the gravitational field around non-Kerr BHs and in which 
the deviations from the Kerr geometry are quantified by 
a set of “deformation parameters”. The simplest non¬ 
rotating BH has only one non-vanishing deformation pa¬ 
rameter and its line element reads 


ds 2 





dt 2 




dr 2 + r 2 dLl , (13) 


where e is the deformation parameter and the gravita¬ 
tional force is weaker (stronger) than the one around 
a Schwarzschild BH with the same mass when e > 0 
(e < 0). If e = 0, we recover the Schwarzschild solution. 
We note that the horizon is at r# = 2 M for any value 
of e. The metric in Eq. (13) is not a solution of any 
known field equations of an alternative theory of gravity. 
It is just the Schwarzschild metric with an ad hoc defor¬ 
mation, which is quantified by the parameter e. Like the 
parameters (3 and 7 in the Parametrized Post-Newtonian 
formalism extensively employed to test general relativity 
in the Solar System, the deformation parameter e in the 
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Johannsen-Psaltis metric is to be thought as a free pa¬ 
rameter to be determined by observations. 

In the case of gravitational wave signals, we can just 
image that the line element in Eq. (131 is the vacuum 
static solution of some alternative theory of gravity. Since 
we want to remain as generic as possible, we do not want 
to specify (and actually we do not know) the exact field 
equations of the underlying theory. We can thus assume 


that the field equations are still given by Eq. ([!]), but 
now is the Johannsen-Psaltis metric. We can then 
proceed as in the previous section and find the funda¬ 
mental equations governing the axial perturbations gen¬ 
erated by a test-particle moving along the geodesic of 
the background metric. Within this approach, axial per¬ 
turbations are still described by Eq. Q, but now the 
potential Vj(r) and the source term Sim\t,r ) are 


w=(i-“ 

Sim(t,r) = - 167 n ( 1 - 


M 3 \ 6M 

1 + e— 5 - )- 5 - 



, 2 ( M 3 \ „ 

1-)(l+ e— 3 - Qlr, 


(14) 

(15) 


Of course, for e = 0 we recover the standard case, namely 
Schwarzschild metric and general relativity. When e 7 ^ 0, 
we have a different wave form Zi m (t) and a different en¬ 
ergy spectrum dW/dtu. In this sense, we have introduced 
a phenomenological parametrization in the general rela¬ 
tivity axial perturbations of a Schwarzschild BH and de¬ 
viations from standard predictions are quantified in terms 
of e. However, such a deformation parameter e is not the 
original e in the Johannsen-Psaltis metric any more. The 
parameter appearing in Eq. (14) and Eq. (15) is a sort of 
effective parameter that takes deviations from the Ein¬ 
stein equations into account. The relation between this 
e and the one appearing in the Johannsen-Psaltis met¬ 
ric and constrained from the electromagnetic spectrum of 
BH candidates [34] is unknown at this level and therefore 
it is difficult to make any comparison between the two ap¬ 
proaches. Of course, we could have also started from the 
very beginning from Eq. Q, writing some ad hoc cor¬ 
rection in terms of a number of deformation parameters. 
However, this might have introduced additional arbitrari¬ 
ness and included unphysical deviations that cannot be 
obtained from any theory of gravity. 


IV. SIMULATIONS 


The numerical method to compute the wave form and 
the energy spectrum is the same as the one adopted in 
Ref. [33]. It consists of three steps. 

1. We fix the parameters of the background metric (e) 
and of the particle ( E and L z ) and we compute the 
particle trajectory for a discrete number of time 
steps, ranging from a large R to the turning point 

R = Rt- 

2 . We compute the source term Si m (t, r) for each time 
and spatial point of the grid. Following Ref. [47], we 
approximate the delta function by a narrow Gaus¬ 


sian 

5{r-R{t))^-°^e - a ^ r - R ^ 2 . ( 16 ) 

V 77 

In our calculations, we use a = 5, which is large 
enough not to affect our results. 

3. We numerically solve Eq. Q. 

With the above machinery, we compute the wave form 
at infinity, Zi m (t), and the energy spectrum, dW/duj , 
for test-particles moving in the gravitational held of 
Schwarzschild and deformed BHs. In our calculations, 
we only compute the signal of the multipole l = 2 and 
m = 1, which is the leading order term for axial pertur¬ 
bations. 

Fig. 0 shows the wave form (top panel) and the en¬ 
ergy spectrum (bottom panel) of the gravitational wave 
signal produced by a test-particle with E = 2.38/r and 
L z = 12.5pM moving in the spacetime of a Schwarzschild 
BH (red solid line) and of deformed non-rotating BHs 
with e = ±2 and 5 (the case e = — 5 has no turning 
point for this choice of E and L z ). To quantify the rel¬ 
ative contribution between the mode excitation and the 
quadrupole orbital emission, we have computed the en¬ 
ergy spectrum with e = 0 in V)(r) and e 7 ^ 0 in S) m (r) 
(quadrupole orbital emission) and then the case with 
e / 0 in b(r) and e = 0 in Si m (r) (mode excitation), 
and compared with the Schwarzcshild and deformed BH 
spectra. The result is shown in Fig. [2j It is evident that 
the gravitational wave signal is mainly determined by the 
geodesic motion of the test-particle in the background 
metric, while the contribution due to the mode excitation 
is smaller. This suggests that in the description of this 
kind of phenomena the field equations of the theory play 
a minor rule with respect to the geodesic motion. This 
may further justify our approach to plug the Johannsen- 
Psaltis metric in the Einstein equations: eventually the 
signal is mainly determined by the quadrupole moment 
formula and the partcle trajectory. The small contribu¬ 
tion from the mode excitation can be easily understood 
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FIG. 1. Wave form (top panel) and energy spectrum (bottom panel) produced by the excitation of axial quasi-normal modes 
of non-rotating BHs (Schwarzschild BH e = 0, deformed BHs e yf 0) by a test particle with E = 2.38 fi and L z = 12.5 fiM. See 
the text for more details. 


if we notice that the maximum of the potential V) (r) is at 
r m ax ~ 3 M and therefore e(M/r max ) 3 -C 1 for e = 0(1). 

Fig-! show the energy spectrum in the case the test- 
particle has different L z . In the top panel we set L z = 
12.0/xM, in the bottom panel L z = 13.0 ^M. We do 
not show the spectrum for e = —2 when L z = 12.0 
because in this case there is no turning point and the 
test-particle is swallowed by the BH. 

In order to figure out whether an accurate measure¬ 
ment of the gravitational wave spectrum can distinguish 
different spacetimes and test general relativity, we can 
proceed as follows. We consider a reference model in 
which the spacetime is described by the Schwarzschild 
metric (e = 0) and the test-particle has a certain spe¬ 
cific energy and axial component of the specific angular 
momentum. Such a reference model can be compared to 


another model, in which the metric has a deformation 
parameter e, and the specific energy and axial compo¬ 
nent of the angular momentum are, respectively, E and 
L z , by evaluating the following function (for the sake of 
simplicity, we set M = fi = 1) 




where ( dW/du ;)* of is the energy spectrum of the reference 
model at the frequency w*, (dW/du) i (e, E, L z ) is the en¬ 
ergy spectrum of the model with parameters (e,E,L z ), 
and C is a co nstant that we have quite arbitrarily set to 
0.3. Eq. ( |7 t| ) clearly looks like a y 2 and C as the error. 
We do not call S y 2 simply because we do not want to 
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FIG. 2. Energy spectrum in the case e = 0 (red solid line), e = 0 in V)(r) but t ^ 0 in Si m (t,r) (green dashed line), t ^ 0 in 
Vz(r) but e = 0 in Si m (t, r ) (blue dotted line), and e ^ 0 (violet dashed-dotted line). Here the test particle has E = 2.38/x and 
L z = 12.5 fiM. See the text for more details. 


perform a rigorous analysis, which would require a more 
detailed discussion beyond the scope of our explorative 
work. 

Fig. a shows the contour level of S assuming that the 
test-particle in the reference model has E = 2.38/r and 
L z = 12.5 M/j,. The sum is performed over 41 frequen¬ 
cies, from w 0 = 0 to W 40 = 0.6/M. In the left panel, S 
is minimized over L z . In the right panel, S is minimized 
over E. If we identified S with y 2 , S = 3.5, 8.0, and 14.2 
would correspond to the probability interval designated 
as 1-, 2-, and 3-standard deviations for three degrees of 
freedom. From this contour levels, it seems that the mea¬ 
surement of the energy spectrum can test the background 
metric and constrain e. There is no correlation between 
the deformation parameter e and the parameters of the 
test-particle E and L z . This is a good result if we want to 


test BH candidates, but we have to consider that we are 
restricting our attention to non-rotating BHs only. The 
typical problem to test BH candidates is the strong corre¬ 
lation between the spin and possible deviations from the 
Kerr solution. We leave the extension to rotating BHs to 
a future work. 


V. SUMMARY AND CONCLUSIONS 

Astrophysical BH candidates are supposed to be the 
Kerr BH of general relativity, but any observational con¬ 
firmation is still lacking and the same Einstein’s theory 
has not yet been tested in strong gravitational fields. 
The Kerr paradigm can be potentially verified by study¬ 
ing the electromagnetic and gravitational signals emitted 
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FIG. 3. As in the bottom panel in Fig. [l] but for a test particle with E = 2.38/r and L z = 12.0/rM (top panel) and with 
E = 2.38 fi, and L z = 13.0/iM (bottom panel). See the text for more details. 


from these systems. The radiation emitted by the gas in 
the inner region of the accretion disk is affected by rela¬ 
tivistic phenomena and the study of the electromagnetic 
spectrum of a BH candidate can thus provide informa¬ 
tion on the background metric close to the compact ob¬ 
ject. Gravitational waves produced by perturbations in 
the spacetime geometry around BH candidates are de¬ 
termined by both the background metric and the field 
equations of the gravity theory. 

In this paper, we have considered the axial per¬ 
turbations generated by a test-particle moving in the 
gravitational field of non-rotating BHs. Starting from 
the Johannsen-Psaltis parametrization, we have derived 
some effective equations governing the evolutions of 
the gravitational wave signal and we have compute 
the associated energy spectrum. Deviations from the 
Schwarzschild predictions in general relativity are quan¬ 


tified by an effective deformation parameter e. Because 
of the dependence of the gravitational wave signal on the 
field equations of the underlying theory, it is not possible 
to directly compare the results from electromagnetic and 
gravitational spectrum. In the former case, experiments 
can really constrain the deformation parameter appear¬ 
ing in the Johannsen-Psaltis metric. With the gravita¬ 
tional wave approach, the final parameter has also ab¬ 
sorbed possible deviations from the Einstein equations. 

We find that the gravitational wave signal is mainly 
determined by the geodesic motion of the particle, while 
deviations in the axial quasi-normal modes due to a dif¬ 
ferent background are smaller. Employing a simple anal¬ 
ysis with the S function in Eq. ( p~T|) , we have compared 
different models to figure out if and how the energy spec¬ 
trum can constrain the deformation parameter e and thus 
test general relativity. Our result is promising: as shown 




















FIG. 4. Contour levels of S. The reference model has e = 0 (Schwarzschild BH), E = 2.38, and L z = 12.5 (fi = A1 = 1). In the 
left panel, S is minimized over L z . In the right panel, S is minimized over E. See the text for more details. 


in Fig. [4j it seems that the measurement of the gravi¬ 
tational wave spectrum can distinguish different space- 
times, constraining e and determining the parameters of 
the test-particle E and L z . However, our discussion is 
limited to non-rotating BHs. The typical problem to test 
BH candidates is the strong correlation between the esti¬ 
mate of the spin and of the deformation parameters. We 
leave the analysis of rotating BHs to a future work, which 
is fundamental to understand if and how this approach 
can test BH candidates. 
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